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ABSTRACT 

We report the discovery of the host galaxy of dark burst GRB 080607 at z G rb = 3.036. GRB 080607 is a 
unique case of a highly extinguished (Ay « 3 mag) afterglow that was yet sufficiently bright for high-quality 
absorption-line spectroscopy. The host galaxy is clearly resolved in deep HST WF3/IR F160W images and 
well detected in the Spitzer IRAC 3.5/.tm and 4.5/im channels, while displaying little/no fluxes in deep optical 
images from Keck and Magellan. The extremely red optical-infrared colors are consistent with the large 
extinction seen in the afterglow light, suggesting that the large amount of dust and gas surface mass density 
seen along the afterglow sightline is not merely local but likely reflects the global dust content across the 
entire host galaxy. Adopting the dust properties and metallicity of the host ISM derived from studies of early- 
time afterglow light and absorption-line spectroscopy, we perform a stellar population synthesis analysis of 
the observed spectral energy distribution to constrain the intrinsic luminosity and stellar population of this 
dark burst host. The host galaxy is best described by an exponentially declining star formation rate of e- 
folding time r = 2 Gyr and an age of ~ 2 Gyr. We also derive an extinction corrected star formation rate of 
SFR ?» 125/i~ 2 M Q yr" 1 and a total stellar mass of M* ~4x 10 11 /i~ 2 M . Our study provides an example of 
massive, dusty star-forming galaxies contributing to the GRB host galaxy population, supporting the notion 
that long-duration GRBs trace the bulk of cosmic star formation. 

Subject headings: gamma rays: bursts individual (080607) — ISM: abundances — ISM: dust, extinction 



1. INTRODUCTION 

Early-time spectra of bright 7-ray burst (GRB) afterglows 
have revealed numerous absorption features that allow accu- 
rate measurements of the chemical composition, dust content, 
and kinematics in the interstellar medium (ISM) of the host 
galaxies (e.g. Fynbo et al. 2006; Savaglio 2006; Prochaska 
et al. 2007; 2008). But as much as 50% of long-duration 
GRBs show a significant suppression in their optical afterglow 
light (Jakobsson et al. 2004; Cenko et al. 2009). While some 
of these "dark" bursts occur during the reionization epoch at 
redshifts z > 6 (e.g. Kawai et al. 2006; Greiner et al. 2009; 
Tanvir et al. 2009; Salvaterra et al. 2009), most result from 
large extinction columns in the ISM surrounding massive star- 
forming regions at more typical redshifts of z = 1 -4 (e.g. Per- 
ley et al. 2009). 

GRB 080607 at redshift zgrb = 3.0363 is a unique case of 
a highly extinguished (Ay ~ 3 mag) afterglow that was yet 
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sufficiently bright for high-quality spectroscopy (Prochaska 
et al. 2009). The afterglow spectrum displays positive detec- 
tions of CO A -X bandheads (Morton & Noreau 1994) that 
have also been seen through translucent molecular gas of the 
Milky Way (e.g. Sonnentrucker et al. 2007). The presence 
of Gell 1602 and OI 1355 absorption features indicates that 
the host ISM has been enriched to roughly solar metallicity. 
Identifications of vibrationally excited H2 indicate the pres- 
ence of substantial molecular gas at a few hundred pc from the 
burst (Sheffer et al. 2009). The large gas surface mass den- 
sity (w 400M Q pc" 2 ) and large molecular gas content found 
in GRB 080607 are unprecedented among either damped Lya 
absorbers along random QSO sightlines or GRB host galax- 
ies (c.f. Srianand etal. 2008; Noterdaeme et al. 2009) Con- 
trary to the common expectation of GRBs occurring prefer- 
entially in low-mass and low-metallicity environments (e.g. 
Fruchter et al. 2006), the observed large metal and dust con- 
tents, together with the mass-metallicity relation known for 
z = 2 - 3 galaxies (e.g. Mannucci et al. 2009), imply that the 
host galaxy is massive and intrinsically luminous. Search- 
ing for the host galaxy of GRB 080607 therefore bears sig- 
nificantly on our general understanding of GRB host galaxies 
and particularly the dark burst population. 

Here we report the discovery of the host galaxy of 
GRB 080607 in an extensive imaging follow-up campaign. 
The observed broad-band spectral energy distribution (SED) 
from optical to IR wavelengths, together with known dust 
properties from studies of early-time afterglow light (Per- 
ley et al. 2010) and known ISM metallicity from absorption- 
line spectroscopy (Prochaska et al. 2009), allow us to con- 
strain the intrinsic luminosity and stellar population of the 
host galaxy of this dark burst. We adopt a A cosmology, 
fiivi = 0.3 and £!a = 0.7, with a dimensionless Hubble constant 
h = //(>/( 100 km s" 1 Mpc" 1 ) throughout the paper. 
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TABLE 1 

Summary of GRB 080607 Host Photometry 



Telescope/Instrument 


Bandpass 


Brightness" 


Keck/LRIS 


G 


AB = 28.2 ±0.4 


Magellan/IMACS 


r 


AB > 27.0 


Keck/LRIS 


I 


AB > 27.0 


HSTAVF3/IR 


F160W 


AB = 24.72 ± 0.03 


Keck/NIRC 


Ks 


AB = 24.8 ±0.7 


Spitzer/IRAC 


3.5 (im 


AB = 20.5 ±0.1 


Spitzer/IRAC 


4.5 /im 


AB = 19.9 ±0.1 


JCMT/SCUBA-2 


850 ^m 


/ < 7.6 mjy 



a When the host is not detected, we measure a 2-<r upper limit for the optical 
and near-IR bandpasses and a 4-<r limit for the 850 fim flux. All magnitudes 
are corrected for Galactic extinction, E(B — V) = 0.023. 



2. MULTI-WAVELENGTH IMAGING OBSERVATIONS AND DATA 
ANALYSIS 

We have carried out an extensive search of the host galaxy 
of GRB 080607 in ground-based and space-based imaging 
observations. Deep optical images of the field around 
GRB 080607 were obtained using the Low Resolution Imag- 
ing Spectrometer (LRIS; Oke et al. 1995) on the 10 m Keck 
I telescope on the night of 2009-Feb-19 UT. Total integra- 
tion times of 2490 s and 2220 s were taken through the G- 
and /-band, respectively. Deep optical r-band images were 
obtained using the short camera in the IMACS multi-object 
imaging spectrograph (Dressier et al. 2006) on the Magellan 
Baade telescope on the night of 2009-May-22 UT. Three ex- 
posures of 600 s each were acquired. Images were flat-fielded, 
registered, and stacked using standard techniques. The final 
stacked images have a mean seeing of 0.8", 0.8", and 1.1" 
in G, r, and /, respectively, and are calibrated using common 
stars in the SDSS photometric catalog. 

Near-IR images of the field were obtained using the Near 
Infrared Camera (NIRC) and the K, filter on Keck I, on the 
night of 2009-May-31 UT A total of 36 exposures of 100 
s each were acquired. Images were processed and stacked 
using standard techniques via a custom Python pipeline. The 
stacked image has a mean seeing of 0.75" and is calibrated 
using standard stars observed throughout the night. 

Deep near-infrared images of the field were obtained us- 
ing the IR channel in Wide Field Camera 3 (WFC3) and 
the F160W filter on board the Hubble Space Telescope (PI: 
Chen). The observations were carried out on 2010-Jul-25 UT. 
Three sets of four exposures (900 s each) were obtained. Indi- 
vidual exposures were reduced using standard pipeline tech- 
niques, corrected for geometric distortion using drizzle, reg- 
istered to a common origin, filtered for deviant pixels, and 
combined to form a final stacked image. 

Infrared images of the field were also obtained on 2010- 
Aug-08 UT by the Spitzer Space Telescope (PI: Perley) in 
both available warm-mission IRAC channels (3.6 and 4.5 
/im). A total of 45 exposures of 100 s each were acquired 
in each channel using a cycling dither pattern. The post- 
BCD calibrated mosaics were retrieved from the Spitzer data 
archive. Finally, we obtained 850 /im images of the GRB field 
under the SCUBA-2 (Holland et al. 2006) Shared Risk Obser- 
vations on the James Clerk Maxwell Telescope (M09BI109, 
PI: Wilson) on 2010-Mar-14 and 2010-Mar-24. A total inte- 
gration time of two hours was acquired. The data were re- 
duced and calibrated as described in Dempsey et al. (2010). 
The combined image has an r.m.s. noise of 2.8 mjy per beam 
with a beam size of 15". 



Optical and IR images of the field surrounding 
GRB 080607 are presented in Figure 1. Astrometric so- 
lutions were obtained using known SDSS objects with a 
mean r.m.s. scatter of « 0.1". The location of the GRB 
afterglow is marked by the cross, which is determined based 
on relative astrometry using seven common stars in our HST 
observations and in early-time afterglow images obtained 
~ 2 hours post burst by UKIRT. The relative astrometry 
provides a precise afterglow position with an error radius 
of 0.05" (c.f. Mangano et al. 2008). At 0.35" away, an 
extended source is clearly detected in the WF3/IR F160W 
image with a half-light radius of r\n « 0.3". We measure 
AB(F160W) = 24.72 ±0.03 over a 0.8"-radius aperture for 
this source. 

To determine whether this source is the host galaxy of 
GRB 080607, we estimate the probability of finding a ran- 
dom foreground galaxy that has its optical disk intercepting 
the afterglow sightline At close projected distances (< 20 h~ l 
kpc), we expect that any foreground galaxy would imprint 
a strong Mg II absorption feature in the afterglow spectrum 
(e.g. Chen et al. 2010). The available afterglow spectrum 
of GRB 080607 allows observations of Mg II absorbers at 
z = 0.895-2.2, and indeed two strong Mgll absorbers have 
been detected at z = 1.341 and z = 1.462 (Prochaska etal. 
2009). At the same time, galaxies A and B are seen within 
3.5" radius of the afterglow sightline with AB(F160W) = 22.4 
and AB(F160W) = 22.0, respectively. Their observed optical 
and near-IR colors are bluer than the host candidate (bottom- 
left panel of Figure 1) and are consistent with galaxies at 
z ~ 1.4. The projected distances are within the expected ex- 
tent of strong Mgll absorbers (Chen et al. 2010). We there- 
fore attribute the two strong Mg II absorbers to galaxies A and 
B. To estimate the probability of a random galaxy at z < 0.895 
occuring within 2 x of the afterglow position, we calcu- 
late the volume density of galaxies with L > 0.025 L* (corre- 
sponding to AB(F160W) = 24.7 at z = 0.895) using the galaxy 
luminosity function of Faber et al. (2007). We find that the 
probability of finding a random z < 0.895 galaxy within this 
small volume is < 1 % which, together with the presence of 
A and B, lead us to conclude that the extended source at the 
afterglow position is the most likely host of GRB 080607. 

The host is also detected in the G-band image with AB(g) = 
28.2 ± 0.4 over a 0.8"-radius aperture, and in the IRAC 
3.5/im and 4.5/im channels with AB(3.5/x) = 20.5 ±0.1 and 
AB(4.5/i) = 19.9 ±0.1 over a 1.2"-radius aperture. The 
adopted aperture sizes roughly matches the size of the aper- 
tures adopted for the optical and near-IR images after account- 
ing for the differences in the PSFs. The errors in the IRAC 
photometry of the host are dominated by contaminating light 
from A and B, and are estimated based on the summed flux 
in a sequence of increasing-diameter apertures. In the /T-band 
image, we observe at the host position a 1.5-u flux detection, 
corresponding to AB(K) = 24.8 ± 0.7, or a 2-cr upper limit of 
AB(K) = 24.5. 

The host galaxy is not detected in the r- and /-band. We 
measure 2-cr upper limits of AB(r) = 21.0 and AB(I) = 27.0. 
The host is not detected in the 850 /im image. Based on 
the deboosted 850 /im fluxes of sub-mm sources observed 
in SCUBA images of comparable r.m.s. noises (Coppin et al. 
2006), we estimate a 4-a upper limit to the 850 /im flux of 7.6 
mjy for the GRB host. 

3. THE HOST GALAXY OF GRB 080607 
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FIG. 1 . — Astrometrically aligned images, 18" on a side of the field surrounding GRB 080607 at zqrb = 3.036. North is up and East to the left. Astrometric 
solutions were obtained using known SDSS objects with a mean r.m.s. scatter of PS 0. 1". The position of the early-time afterglow is at the cross with 0.05" 
error. The ground-based images have been smoothed over the FWHM of the PSF for revealing faint features. An extended source is clearly detected at 0.35" 
away from the afterglow position in the HST WF3/IR F160W image (arrow in the top right panel), which is also well detected in the Spitzer IRAC 3.5/xm and 
4.5^tm images (two bottom right panels). While the souce is not seen in the ground-based r and / images, it displays traces of flux in the deep G- and S'-band 
images. The consistent red color (bottom-left panel) and astrometric match with the afterglow lead us to conclude that this extended source is the most likely host 
galaxy, H 080607. Sources A and B exhibit bluer colors and are likely the absorbing galaxies of two strong Mgll absorbers (Prochaska et al. 2009) at z = 1.462 
and z = 1.361, respectively. 



To constrain the stellar population and star formation his- 
tory of the host galaxy, we consider a suite of synthetic stel- 
lar population models generated using a revised Bruzual & 
Chariot (2003) spectral library that includes a new prescrip- 
tion for the TP-AGB evolution of low- and intermediate-mass 
stars (Marigo & Girardi 2007). To improve the uncertainties 
in the model analysis (such as the age-metallicity degeneracy, 
e.g. Worthey 1994), we take into account known dust prop- 
erties from studies of early-time afterglow light (Perley et al. 
2010) and known ISM metallicity from absorption-line spec- 
troscopy (Prochaska et al. 2009). 

The SED of the host galaxy is found to be extremely red, 
consistent with the large extinction seen in the afterglow light. 
As demonstrated in Perley et al. (2010), the observed SED of 
the GRB afterglow is best-described by a combination of a 
single power-law spectrum (characteristic of the intrinsic af- 
terglow radiation) and an extinction law that is characterized 
by the broad 2175-A absorption band, commonly seen in the 
Milky Way (MW) and the Large Magellanic Cloud (LMC), 
and a large extinction column of Ay ~ 3.1 at z = 3. How- 
ever, the best-fit extinction curve appears to be relatively flat 
at UV wavelengths with Ry ~ 4, and the 2175-A absorption 
feature is not as strong as what is seen in MW or LMC. Conse- 
quently, Perley et al. (2010) derived a best-fit extinction curve 
for the host using a general extinction law from Fitzpatrick & 
Massa (1990), which allows variations in the UV slopes, and 
the depth and width of the 2175-A absorption. 

The large Ry suggests that the GRB occurred in a dense 
environment (e.g. Valencic et al. 2004), consistent with the 
large gas column seen in afterglow spectra. The extremely 
red colors of the host galaxy suggest that the large amount of 
dust and gas surface mass density seen in the afterglow light 
is not merely local to the line of sight to the burst, but likely 
reflects the global dust content across the entire host galaxy. 
In addition, the presence of the 2175-A dust feature suggests 
that the host galaxy resembles mature galaxies like MW or 
LMC, rather than young star-forming systems like the Small 
Magellanic Cloud. In the following stellar population syn- 
thesis analysis, we include priors from the best-fit Fitzpatrick 
& Massa extinction law {Ry = 4. and Ay = 3) of Perley et al. 
(2010) and solar metallicity measured for the host ISM from 
afterglow absorption-line observations. We also examine pos- 
sible biases due to adopting these priors. 

For the stellar population library, we adopt a Chabrier initial 
mass function (IMF) with a range of star formation histories 
from a single burst model to an exponentially declined star 



formation rate (SFR) model with an e-folding time ranging 
from t = 0. 1 - 2 Gyr. Over the spectral range from rest-frame 
~ 1000 A to ~ 1 micron, the observed emission is dominated 
by stellar light in dusty galaxies (Hainline et al. 2009). We 
therefore consider only stellar components in the SED mod- 
els. We perform a maximum likelihood analysis to compare 
the observed SED and a grid of model expectations, taking 
into account both detections and non-detections. The likeli- 
hood function of this analysis is defined as 

r, ^ (ri S 1 age) ] 2 1^ 
£(r,age)= M |exp<^ -- Hx 

where fj is the observed flux (in fiJy) of the host in band- 
pass i, f is the model expectation, and cr, is the measurement 
uncertainty of fi. The first product of Equation (1) extends 
over the k measurements and the second product extends over 
the / upper limits. The result shows that, with the adopted 
solar metallicity and a relatively grey extinction law, the stel- 
lar population of the host galaxy is best described by an ex- 
ponentially declining SFR of r = 2 at the age of ~ 2.2 Gyr. 
The observed SED and the best-fit model are presented in the 
top panel of Figure 2, showing a good agreement in all band- 
passes. The likelihood function of the stellar age is presented 
in the bottom panel of Figure 2. 

The best-fit stellar age of the host is similar to the age of the 
universe at z = 3, indicating that the host galaxy was formed 
very early in time. The comparable stellar age and star for- 
mation e-folding time indicates that the host has been un- 
dergoing active star formation since birth. We estimate the 
SFR and total stellar mass of the host galaxy based on ex- 
tinction corrected UV luminosity and mass-to-B-band light 
ratio. The unobscured best-fit spectrum is displayed in the top 
panel of Figure 2 (dash-dotted curve). Adopting the best-fit 
model, we derive a rest-frame, extinction corrected UV lumi- 
nosity of L(1500)= 1.8 x lO^/z^ergs" 1 Hz" 1 at 1500 A. For 
a Chabrier IMF, this corresponds to SFR « 125/r 2 M Q yr -1 
(Salim et al. 2007). At the age of « 2 Gyr, we derive a total 
stellar mass of M* ~ 4 x 10 11 h~ 2 M . The uncertainty in M* 
due to uncertainties in the star formation history is found to 
be w 50%. 

To investigate possible biases due to adopting the afterglow 
extinction curve and line-of-sight metallicity for the global 
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FIG. 2. — Top: Comparison of the observed SED of the dusty host of 
GRB 080607 at zgrb = 3.036 and the best-fit stellar population synthesis 
model. Optical and infrared photometric measurements are shown in open 
circles with errorbars. We place a 2-cr upper limit to the observed bright- 
ness, when the galaxy is not detected. The horizontal errorbars indicate the 
FWHM of each bandpass. The solid curve represents the best-fit synthetic 
model after accounting for dust obscuration (the best-fit Fitzpatrick & Massa 
law of Perley et al. 2010 with Ry = 4 and Ay = 3 at z = 3). The open squares 
represent the predicted brightness from the best-fit model. The dash-dotted 
spectrum at the top shows the intrinsic spectrum prior to the application of 
dust obscuration. The dotted curve represents a best-fit model based on the 
same extinction law but with Ay =2.5, indicating that reducing the amount of 
dust extinction no longer provides a good fit to the observed SED. Adopting 
the 2-(T upper limit for the S'-band, AB(K) > 24.5, does not change the re- 
sults. Bottom: The likelihood function of the minimum age of the underlying 
stellar population as described by the broad-band SED. 

properties of the host galaxy, we repeat the stellar population 
synthesis analysis with varying metallicity and the amount of 
dust extinction. We find that for a fixed amount of extinc- 
tion, the derived SFR and are insensitive to the adopted 
metallicity. Reducing the amount of dust extinction in the 
model SED results in a poor fit to the observed SED for any 
combination of star formation history and metallicity (e.g. the 
dotted curve in Figure 2). This exercise confirms that the dust 
extinction law determined from the afterglow light and the 
line-of-sight metallicity are representative of the mean prop- 
erties across the entire host galaxy. We therefore conclude 
that the best-fit SFR and are robust. 

4. DISCUSSION 

Our multi-wavelength imaging follow-up has uncovered 
an extremely red galaxy at the location of the "dark" 
GRB 080607 at zgrb = 3.036. Given the coincident position 



and consistent red optical and near-infrared colors, we argue 
that the galaxy is the host of the dusty burst. The host galaxy 
is clearly extended in the HST WF3/IR F160W image with a 
half-light radius of fs 0.3", corresponding to a physical half- 
light radius of 1.8 h~ l kpc at z = 3. This is typical of what 
is seen for star-forming galaxies at z ~ 3 (e.g. Bouwens et al. 
2004). The host galaxy is both massive and actively form- 
ing stars. We estimate the mean radiation field Iq in the host 
ISM at near-UV wavelengths (w 1500-2000 A) by averag- 
ing the extinction-corrected UV flux over the extent of the 
host galaxy seen in the HST WF3/IR image, and find that 
/o«3x 10~ 4 photons cm" 2 s" 2 Hz" 2 . The estimated mass and 
SFR are among the highest known in the GRB host galaxy 
population (e.g. Christensen et al. 2004; Savaglio et al. 2009; 
Chen et al. 2009) and in star-forming galaxies at z ~ 3 (e.g. 
Erb et al. 2006), but are comparable to the average of sub- 
mm sources (e.g. Borys et al. 2005; Dye et al. 2006; Serjeant 
et al. 2008) and more than five times lower than the brightest 
sub-mm galaxies (e.g. Tacconi et al. 2006). 

Unobscured GRB host galaxies at intermediate redshifts ap- 
pear to be underluminous and low mass systems (e.g. Fruchter 
et al. 2006). It has therefore been argued that long-duration 
GRBs originate preferentially in relatively metal deficient 
star-forming regions (Wolf & Podsiadlowski 2007; Modjaz 
et al. 2008). A low-metallicity environment is favored by pop- 
ular progenitor models so that the progenitor star can preserve 
high spin and a massive stellar core to produce a GRB (e.g. 
Yoon & Langer 2005; Woosley & Heger 2006). However, 
statistical studies have shown that the distributions of ISM 
metallicity and UV luminosity of known GRB host galaxies at 
z > 2 are consistent with the expectations of a UV luminosity 
selected star-forming galaxy sample (Fynbo et al. 2008; Chen 
et al. 2009). 

GRB080607 represents a unique example of a dark GRB 
that was luminous enough to allow detailed observations of 
its afterglow despite it occurring in a heavily obscured galaxy. 
Thus we were able to secure an unambiguous redshift of the 
dusty burst through afterglow absorption-line spectroscopy 
and to identify the host galaxy based on an astrometric match 
to the afterglow position. Our current limit at 850/im is typ- 
ical of what is seen in GRB host galaxies (e.g. Tanvir et al. 
2004). Improved sensitivity limit at the sub-mm wavelengths 
will provide a different estimate of the SFR (e.g. Chapman 
et al. 2005) and some constraint on the dust-to-stellar mass 
ratio of the host. The large gas and dust mass uncovered in 
the afterglow spectrum together with the IR bright host galaxy 
already show that mature, dusty star-forming galaxies do con- 
tribute to the GRB host galaxy population (see also Levesque 
et al. 2010), supporting the notion that long-duration GRBs 
trace the bulk of cosmic star formation. Follow-up near-IR 
spectroscopy of the host will not only confirm the host identi- 
fication, but also allow a detailed study of the gas kinematics 
in the host ISM. 

We thank John Mulchaey for obtaining the IMACS r im- 
ages presented in this paper and Stephane Chariot for pro- 
viding the updated spectral library. Support for the HST DD 
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Erratum: "A Mature Dusty Star-forming Galaxy Hosting GRB 080607 at z = 3.036" (2010, ApJ, 723, L218) 

We have discovered an error in the photometric measurements of the host galaxy in our Spitzer IRAC images. The host is 
detected in the IRAC 3.5^m and 4.5/mi channels with Afi(3.5/x) = 22.9 ± 0.2 and AB(4.5^) = 22.7 ± 0.2 mag. The photometric 
measurements of the host galaxy in other bandpasses remain unchanged. Adopting the revised Spitzer IRAC photometry and the 
original optical and near-IR photometric measurements, we estimate the total stellar mass (M„) and on-going star formation rate 
(SFR) of the host galaxy based on the stellar population synthesis analysis described in Chen et al. (2010). Given the uncertainties 
in the global dust content of the host galaxy, we allow the metallicity, Ay, and dust extinction law to differ from what were found 
in the afterglow light (e.g. Prochaska et al. 2009; Perley et al. 2010). The likelihood analysis described in Equation (1) of Chen 
et al. (2010) yields an extinction-corrected SFR of (8-12) h~ 2 M Q yr" 1 , a mean ISM radiation field 7 « (2.3-3.5) x 10" 5 photons 
cm" 2 s" 2 Hz" 1 , and M* = (0.5- 1.4) x 10 10 /r 2 M Q for the host galaxy. These are about an order of magnitude lower than those 
originally published. We note that the uncertainties in the derived M», Iq and SFR are driven by the uncertainties in the global 
dust extinction law of the host galaxy. The galaxy is still fairly massive, but not as extreme as previously thought. The observed 
spectral energy distribution (SED) of the host galaxy is presented in the revised figure below, together with the best-fit synthetic 
model of super-solar metallicity and a Milky-Way type dust extinction law of Ay = 1.25. A Fitzpatrick & Massa (FM) law 
described in Perley et al. (2010) with Ay = 1.8 and super-solar metallicity produces a similarly good fit to the observed SED. 

Adopting the priors from the best-fit dust obscuration, the Fitzpatrick & Massa (FM) law from Perley et al. (2010) with Ry = 4.2 
and Ay = 3.3 at z = 3 and solar metallicity from Prochaska et al. (2009) leads to a similar estimate of M» (~ 5 x 10 9 /i~ 2 M Q ) but 
significantly higher SFR (~ 230 /T 2 M Q yr" 1 ). However, this model also predicts an observed optical brightness that is AAB = 1 
mag brighter than the observed 2-a upper limits in the r and / bands. We therefore consider this model unlikely to represent 
the global extinction property of the host galaxy. The result indicates a large spatial variation in the dust content across the host 
galaxy. 




FIG. 3. — Top: Comparison of the observed SED of the dusty host of GRB 080607 at zgrb = 3.036 and the best-fit stellar population synthesis models. Optical 
and infrared photometric measurements are shown in open circles with errorbars. We place a 2-cr upper limit to the observed brightness, when the galaxy is not 
detected. The horizontal errorbars indicate the FWHM of each bandpass. The solid curve represents the best-fit synthetic model of super-solar metallicity and a 
Milky-Way type dust extinction law of Ay = 1.25. The open squares represent the predicted brightness from this best-fit model. The thin dash-dotted spectrum at 
the top shows the intrinsic spectrum prior to the application of dust obscuration. Bottom: The likelihood functions of the minimum age of the underlying stellar 
population as described by the broad-band SED. 



